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Abstract 

We present a scenario for building the equatorial ridge and de-spinning lapetus through 
an impact-generated disk and satellite. This impact puts debris into orbit, forming a ring 
inside the Roche limit and a satellite outside. This satellite rapidly pushes the ring material 
down to the surface of lapetus, and then itself tidally evolves outward, thereby helping to 
de-spin lapetus. This scenario can de-spin lapetus an order of magnitude faster than when 
tides due to Saturn act alone, almost independently of its interior geophysical evolution. 
Eventually, the satellite is stripped from its orbit by Saturn. The range of satellite and 
impactor masses required is compatible with the estimated impact history of lapetus. 



1. Introduction 



The surface and shape of lapetus (with equatorial ra dius, -Rt=746 krn, and bulk density, 
p = 1.09 g cm~^) are unlike those of any other icy moon ( Jacobson et al. 2006h. About half 



see 



Porco et al. 2005. for 



of lapetus' ancient surface is dark, and the other half is bright 
discussion). This asymmetry has been explained recently as the mi gration of water ice due t o 
the deposition of darker material on the leading side of the body ( Spencer and Den^ 201oh. 



lapet us also has a ridge system near its equator, extendin g > 110° in longit ude (jPorco et al. 
20051 ) . that rises to heights of ~ 13 km in some locations (IGiese et a/.ll2008l ). The ridge itself 



is heavily cratered, suggesting it originated during lapetus' early history. Finally, lapetus' 
present-day overall shape i s consistent with a rapid 16-hour spin period rather t han its 
present 79-day spin period (lThomasll2010l : ICastillo-Rogez et a/.l 120071 : lThomasll2010l ). 



To some, the equatorial position of the ridge and lapetus' odd shape suggest a causal 
relationship. Most current explanations invoke endogenic processes. For example, detailed 
models of lapetus' early thermal evolution suggest that an early epoch of heating due to 
short-lived ^^Al and ^°Fe is required to close off primordial porosity in the object while 
simultaneously allowing it to rapidly de-spin, cool, and l ock in a "fossil bulge " indicative 
of an early faster spin period (jCastillo-Rogez et a/.l 120071 : iRobuchon et a/.ll2010l ). Recently, 
Sandwell & Schubert (2010) suggested a new and innovative mechanism for forming the 
bulge and ridge of lapetus through contraction of primordial porosity and a thinned equa- 
torial lithosphere. However, only a narrow range of parameters allows both a thick enough 
lithosphere on lapetus to support the fossil bulge, while also being sufficiently dissipative to 
allow lapetus to de-spin due to Saturn's influence on solar system timescales. 

In these scenarios, the ridge represents a large thrust fault arising from de-spinning. One 
difficulty faced by these ideas is that the stresses arising from de-spinning at the equato r 
are perpendicular to the orientation required to creat e an equatorial ridge (IMeloshl Il977l). 
Othe r interior processes, such as a convective upwelling (jCzechowski and Leliwa-Kopvstvhskil 
2008 ). or convection coupled with tidal dissipation driven by the de-spinning (IRoberts and Nimmo 
20091 ) are required to focus and reorient de-spinning stresses on the equator. These latter 



- 3 - 



models have difficulty reproducing the ridge topography be cause thermal buoyancy str esses 
are insufficient to push the ridge to its observed height (see lDombard and Cheng |2008| ). 



Alternatively, the ridge may be exogenic. One leading hypo t hesis is that the ridg e 
represents a ring system deposited onto lapetus' surface (llpl l2006l : iDombard et a/.l 120101 ) . 
This model has the benefit of providing a natural explanation for the mass, orientation, and 
continuity of the ridge, which present a challenge to endogenic models. 

Here we extend this idea to include a satellite that accretes out of the ring system 
beyond the Roche limit. As we show below, this can significantly aid in the de-spinning of 
lapetus. In particular, we hypothesize that: 

1) lapetus suffered a large imp act that produced a debris disk similar to what is believed 
to have formed Earth's Moon JCanupl[20M : llda et almWA iKokubo et o/lboooh . Like the 
proto-lunar disk, this disk straddled the Roche radius of lapetus, and was quickly coUisionally 
damped into a disk. As a result, a satellite accreted beyond the Roche radius, while a 
particulate disk remained on the inside. Also, the impact left lapetus spinning with a period 
< 16 hr, thereby causing the bulge to form0 

2) Gravitational interactions between the disk and lapetus' satellite (hereafter known as the 
sub-satellite) pushed the disk onto lapetus' surface, forming the ridge. As Ip (2006) ffist 
suggested, a coUisionally damped disk, similar to Saturn's rings, will produce a linear feature 
precisely located along the equator. Thus, it naturally explains the most puzzling properties 
of the ridge system. The impact velocity of the disk particles would have been only ~ 300 m 
s^^ and mainly tangential to the surface, so it is reasonable to assume that they would not 
have formed craters, but instead piled up on the surface. 

3) Tidal interactions between lapetus and the sub-satellite led to the de-spinning of lapetus 
as the sub-satellite's orbit expanded. Eventually, the sub-satellite evolved far enough from 
lapetus that Saturn stripped it away. lapetus was partially de-spun and continued de- 
spinning under the influence of Saturn. Finally, the sub-satellite was either accreted by one 
of Saturn's regular satellites (most likely lapetus itself) or was ejected to heliocentric orbit 
(cf. The end-state is a de-spun lapetus that has both a bulge and a ridge. Faster 
de-spinning aided by the presence of a sub-satellite likely relaxes constraints on the early 
thermal evolution of lapetus determined by prior works (Castillo- Rogez et al 2010; Robuchon 
et al. 2010). 

Because the results of one part of our story can be required by other parts, we begin our 
discussion in the middle and ffist find, through numerical simulations, the critical distance 



^It is important to note that the impact that we envision is in a region of parameter space that has yet 
to be studied. Such an investigation requires sophisticated hydrodynamic simulations and thus is beyond 
the scope of this paper. We leave it for future work. We emphasize, however, th at the genera l geometry we 
envision has been seen in many hydrodynamic simulations of giant impacts (e.g.. ICanudl2004l) . so we believe 
that our assumed initial configuration is reasonable. 
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(ttst) at which a sub-satelhte of lapetus will get stripped by Saturn. Knowing this distance, 
we integrate the equations governing the tidal interactions between both Saturn and lapetus, 
and between lapetus and the sub-satellite, to estimate limits on the mass of the sub-satellite. 
We then study the fate of the sub-satellite once it was stripped away from lapetus by Saturn. 
Finally, using crater scaling relations we reconcile a sub-satellite impact with the topography 
of lapetus. 



2. Satellites Stripped by Saturn 

The distance at which a satellite of lapetus becomes unstable is important for calculating 
tidal evolution timescales. In systems containing the Sun, a planet, and a satellite, prograde 
satellites are not expected to be stable beyond ~ -Rh/2, where the Hill radius is defined as 
Rn = a(m/3 MY^^ with a as the planet's semi-major axis, m as its rnass, and M as the tota l 



system mass (IHamilton and BurnsI Il99ll : iBarnes and O'BrienI 120021 : iNesvorny et al\ |2003| ) . 



In our case, lapetus plays the role of the planet, and Saturn the role of the Sun. However, 
the tidal evolution timescale depends strongly on semi-major axis (as the —13/2 power, 
Eq. [3]) and thus the success of our model depends sensitively on the value of the critical 
distance, Therefore, we performed a series of numerical simulations to determine agt- 

This experiment used the swift_WHM integrator (Levison & Duncan 1994; which is 
based on Wisdom & Holman 1991) to integrate two sets of test particles consisting of 500 
objects, each of which were initially on orbits about lapetus with semi-major axes, a, that 
ranged from 0.1-0.8 -Rh- The particles in the first set were initially on circular orbits in the 
plane of lapetus's equator. Particles in the second set had initial eccentricities, e, of 0.1, and 
inclinations, i, that were uniformly distributed in cos (i) between i = and i = 15°. Saturn is 
by far the strongest perturber to the lapetus-centered Kepler orbits and is the main source 
of the stripping. For completeness, we have also included the Sun and Titan. The effects 
of the other Saturnian satellites are at least two orders of magnitude smaller than those of 
Titan and thus can be ignored. 

The simulations were performed in an lapetus-centered frame. The lifetime of particles 
dropped precipitously beyond 0.4 _Rh, suggesting that any sub-satellite with a larger semi- 
major axis would very quickly go into orbit around Saturn (Fig. [1]). Thus, we adopt this 
limit, which is equivalent to 21 Ri, in our calculations below. 

[Fig. [H here] 



3. Tidal evolution of lapetus 



The de-spinning of lapetus by Saturn has long been considered problematic, because 
for nominal Q/\k2\ (~10^), lapetus should not have de-spun over the age of the solar system 
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( jPeald Il977l ). Starting with the assumption of constant Q/|A;2|, and using the standard 



de-spinning timescale (Murray & Dermott 1999, eq 4.163), 



"I = ~sigii(iq — n) — — - — I — \ n 



2aQ mi(m|^ + mi) 



a 



where a < 2/5 is the moment of inertia constant of lapetus, mi is its mass, fli is its spin 
frequency, \k2\ is the magnitude of the ^2 Love number, Q the tidal dissipation factor, m^ is 
the mass of Saturn, and a and n are the semi-major axis and mean motion of lapetus. The 
\k2\ and Q values used throughout are for lapetus only. For the tidal interaction between 
lapetus and Saturn, Qi > n, so the effect is always to decrease the spin of lapetus. 

For these simple assumptions, the de-spinning from 16 h to a rate synchronous with the 
orbital period, 79.3 days, takes 3.6 x 10^ (Q/|^2|) years, nominally 36 Gyr, for a density p = 
1 g cm"^. Using detailed geophysical models, Castillo- Rogez et al. (2007) and Robuchon et 
al. (2010) showed Saturn can de-spin lapetus on solar system timescales, although only for 
a narrow range of thermal histories. Our goal here is to investigate how the addition of the 
sub-satellite affects the de-spinning times. 

Given that detailed models of Castillo-Rogez et al. (2007) and Robuchon et al. (2010) 
used different methods, and that we are only interested in how the de-spinning timescale 
changes with the addition of a satellite, we take a simple approach of integrating a modified 
version of Eq. [H Our first adjustment is to remove the assumption of constant Q/\k2\. This 
ratio is dependent on the tidal frequency, {fl — n), and accounts for the manner in which a 
material or body reacts to tidal stresses. We start with a model of lapetus consisting of a 
time-invariant 200-km thick lithosphere with a Maxwell viscoelastic rheology with rigidity 
/i = 3.6 X 10^ P a and viscosity y = 10^^ Pa s , which is strong enough to support the equatorial 



bulge and ridge ICastillo- Rogez et al\ ( 120071 ). overlying a mixed ice/rock mantle with a lower 
viscosity, representative of an interior warmed by radiogenic heating. We performed two 
types of simulations. In the first, the viscosity of the mantle is held constant with time 
and has values from rj = 10^^- 10^*^ Pa s (typical for the interior of an icy satellite at 240 
- 270 K). In the second, we allow t] of the inner ice/rock mantle to vary according to 
the thermal evolution models in Gastillo- Rogez et al. (2007) and Robuchon et al. (2010). 
In particular, we employ the LLRI model of Castillo-Rogez et al. (2007), and the 0.04 and 
72 ppb ^^Al cases from Robuchon et al. (2010). Love numbers are calculated for a spherically 
symmetric, uniform-density lapetus using the SatStress software package (Wahr et al. 2009). 
We calculate the Love number k2 (which is a complex number for a viscoelastic body, see 
Wahr et al., 2009 for discussion) and estimate Q/\k2\ = l/Im(A;2) (Segatz et al. 1988). The 
values of Q/\k2\ vary over an order of magnitude for each value of rj for the important range 
of tidal frequencies. 

An integration of equation ([1]) was performed using a Bulirsch-Stoer integrator for times 
up to 100 Gyr, incorporating the frequency dependent Q/\k2 \ for different internal viscosities 
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which, in turn, is a function of temperature. Without the sub-satelhte, the time for lapetus 
to reach synchronous rotation ranged between 5 xlO^ (fixed r] = 10^^ Pa s) to 2 xlO^^ years 
(0.04 ppb ^^Al case from Robuchon et al. 2010). We describe an investigation of the effect 
that a sub-satelhte could have on the spin of lapetus in the next subsection. 



3.1. Tidal interaction with a sub-satellite 

The sub-satellite raises a tidal bulge on lapetus, causing lapetus to de-spin and the 
sub-satellite's orbit to change. The change in spin rate for lapetus due to a sub-satellite is 
(Murray & Dermott 1999, eq. 4.161), 

3\k2\ ml fRi^^ 



r]i = -sign(a-n) ;''^;' — -(^) (2) 

and, the change in the satellite's orbit by (Murray & Dermott 1999, eq. 4.162), 



3|/i;2| "^ss / -Ri 
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d = sign(ni - n)— -— — na . (3) 

2aQ mi \ a J 

Together Eqs. [2] and [3] describe the interaction between the sub-satellite and lapetus, where 
m^s is the mass of the sub-satellite. The term sign(f2i— n) is of great importance, determining 
whether the satellite evolves outward while decreasing the spin of lapetus, or inwards while 
increasing the spin of lapetus. At semi-major axis Osync = (^(mi + m^s) /^lY^"^ ■, fii = n, 
representing a synchronous state. If the sub-satellite has a < ctsync; it evolves inwards; if 
a > cisync, it evolves outwards. Saturn is gradually decreasing the rotation rate of lapetus, 
and thus the synchronous limit slowly grows larger, possibly catching and overtaking a slowly 
evolving sub-satellite. The orbital period of a sub-satellite at 21 Ri, the distance at which 
we consider a satellite stripped by Saturn, is ~ 12.8 days. Thus, if lapetus is de-spun to a 
period of 12.8 days before the sub-satellite reaches 21 Ri, it will be caught by the expanding 
synchronous limit. 

For the integrations of the sub-satellite's tidal evolution, the sub-satellite's mass is used 
as a free parameter, while the starting semi-major axis is set to 3 Ri. This distance is derived 
from the expected origin of the sub-satellite accretin g from an impact-caused debris disk 



encircling lapetus ( llda et a/.lll997l : iKokubo et a/.ll2000l ). However, tidal evolution timescales 



for a sub-satellite are largely insensitive to the initial semi-major axis, so this starting point 
only needs to be beyond the synchronous limit for the model to be accurate. With a rotation 
period of 16 h for lapetus, agync would have been 2.94 Ri, which is outside the Roche hmit 
defined to be at rroche = 2.46i?i(pi/pss)*-^^^'* ~ 2.53i?i for = 1 g cm~^. Thus, a satellite 
forming at 3.0 Ri would be above the synchronous limit and destined, initially, to evolve 
outward due to tidal interaction with lapetus. 
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Equations (2) and (3) were then integrated with equation (1), to follow the evolution of 
the spin of lapetus due to Saturn and the sub-satellite. We studied the geophysical models for 
lapetus described in the last section, along with sub-satellites with mass ratios, q = rngs/mi, 
between 0.0001 and 0.04. We summarize the suite of simulations in Fig. 2, showing the 
time of de-spinning for lapetus and the time at which the sub-satellite is stripped by Saturn 
or tidally evolves back to re-impact lapetus. We present the data scaled to the de-spinning 
time due to Saturn alone to highlight the effect of the sub-satellite in accelerating the tidal 
evolution of the system. The fate of the system, in regards to the escape or re-impact of the 
sub-satellite and the de-spinning time of lapetus, are separated into three distinct classes of 
outcomes based on q. 

3.1.1. Synchronous lock and re-impact: q > 0.021 

Above a mass ratio q > 0.021, the sub-satellite does not reach Ogt before becoming 
synchronous with the spin of lapetus (see Fig. With both the sub-satellite and Saturn 
working to slow the spin of lapetus, the synchronous limit grows to 21 /2i before the sub- 
satellite evolves to that semi-major axis. This result only varies mildly for the different 
geophysical models of lapetus, as all three timescales depend linearly on Q/\k2\] therefore, 
the re-impact outcome only depends on mass ratio. However, the time to reach this outcome 
ranges from 10^^ yr for the Robuchon et al. (2010) models to 10^ yr for the Castillo-Rogez 
et al. (2007) LLRI model. 

Upon reaching synchronous rotation, the evolution does not stop because Saturn is 
still tidally interacting with lapetus. As Saturn continues to slow the spin rate of lapetus, 
the synchronous limit moves beyond the sub-satellite, which then begins to tidally evolve 
inwards. The sub-satellite is doomed to evolve inwards and hit lapetus. Given that the 
sub-satellite started at 3 Ri and finishes by impacting lapetus, it makes a net contribution 
to the angular momentum of lapetus, and so it is spinning faster than if the sub-satellite 
were never there. Thus, the de-spinning of lapetus (after re-impact) finishes later than it 
would have by Saturn tides alone (see Fig. [3^). 

3.1.2. Satellite is stripped: 0.006 <q< 0.021 

For 0.006 < q < 0.021, the sub-satellite evolves to 21 Ri and is stripped by Saturn 
before attaining a synchronous orbit. As it moves out, the sub-satellite carries away angular 
momentum from lapetus, allowing it to rapidly de-spin. This angular momentum is then 
removed from the lapetus system when the sub-satellite is stripped. 

The sample evolution for a system with a constant rj = 10^^ Pa s and q = 0.018 (Fig.[3]D) 
shows that the sub-satellite reaches an orbital period of ~ 12 days (a = 21 Ri) before lapetus 
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reaches that spin period. It is important to note that the sub-satelhte does not totally de- 
spin lapetus. Instead, it slows lapetus down enough that Saturnian tides (which are faster 
because Q/\k2 \ is a decreasing function of the spin rate) can finish the job. For all but one (see 
§3.1.4p of our geophysical models, the sub-satellite can de-spin lapetus an order of magnitude 
faster than it is when de-spun by Saturn alone. This order of magnitude difference means 
that lapetus could have de-spun in 500 Myr, for situations that would otherwise require the 
age of the solar system. 

3.1.3. Slow evolution of a small satellite: q < 0.006 

The tidal evolution timescale of the sub-satellite's orbit expansion depends on g, and for 
smaller mass ratios, the evolution takes longer. Below q < 0.006, the de-spinning of lapetus 
due to Saturn is fast enough that the location of synchronous rotation sweeps past the sub- 
satellite (see Fig. After this occurs, the sub-satellite is then below the synchronous 
limit and doomed to evolve back in towards lapetus. In this scenario, the evolution of the 
sub-satellite back to the surface of lapetus takes longer than it takes for lapetus to de-spin. 

lapetus de-spins faster than Saturn otherwise could do alone. In this case, however, the 
relevant time constraint becomes the sub-satellite impact time (the small open symbols in 
Fig. rather than the de-spinning time because we currently do not see a satellite. We 
find that sub-satellite impact time can be shorter than the de-spinning time for q > 0.003. 
However, it should be noted that this dynamical pathway only helps by at most roughly a 
factor of 2 over Saturn acting alone. In addition, the sub-satellite is likely to tidally disrupt 
on its way in, forming a second, significantly fresher, ridge. This is probably inconsistent 
with the ridge's ancient appearance. Thus, we think that this particular dynamical pathway 
can probably be ruled out, but we include it for completeness. 

3.1.4. Discussion of sub-satellite tides 

The integrations have bracketed the possible behavior of the Saturn-Iapetus-sub-satellite 
system. At high and low mass ratios the sub-satellite is doomed to return to re-impact 
lapetus, while for 0.006 < q < 0.021 the sub-satellite is stripped. As Fig. [2] shows, sub- 
satellites with masses between 0.005 < q < 0.021 decrease the despinning time over that of 
Saturn alone. This effect can be as large as a factor of 10 for q ~ 0.02. Consistent with prior 
work, we find that a low-viscosity interior (presumably warmed by radiogenic heating) is 
required to despin lapetus over solar system history. We find that the age of the despinning 
event with and without a sub-satellite are similar if the thermal evolution of lapetus follows 
the LLRI model of Castillo- Rogez et al. (2007). In this case, lapetus' interior heats slowly 
during the first Gyr of its history. When the interior is warmed close to the melting point, 
tides drive rapid despinning. The presence of the sub-satellite can significantly shorten 
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this period of time, but because the de-spinning time is short irrespective of whether the 
sub-satelhte is present, it does not significantly alter the time in solar system history when 
lapetus is de-spun. 

4. The ridge and sub-satellite formation 

Given the sub-satellite masses which can assist in de-spinning lapetus, and estimates on 
the mass of the equatorial ridge, a constraint can be placed on the amount of mass placed 
into orbit by the original disk-forming impact. The ring of debris which collapses to form 
the ridge must do so rapidly, as a ring is not currently observed, and the ridge is one of the 
oldest features on the surface of lapetus (Giese et al. 2008). 

4.1. Ridge mass 

The ridge has an unknown mass due to incomplete imaging and significant damage 
from cratering. Ip (2006) estimated its mass assuming that it had, at one time, completely 
encircled the equator with a height of 20 km and width of 50 km, mridgc = 2 x tt x /?i x 
20 km X 50 km with a density of 1 g cm^'^ equals a mass of 4.4x10^^ grams (Ip 2006 used 
a radius of 713 km for this mass estimate). Giese et al. (2008) found a maximum height of 
13 km in Digital Terrain Models (DTM) models, though the shape model maximum height 
was 20 km. Comparing the dimensions given by Ip (2006) with the profiles in Giese et al. 
(2008), we set a lower limit by taking a factor of two in both vertical and horizontal extent 
and assuming that the cross section is a triangle rather than a rectangle, yielding an estimate 
~8 times lower, 5.5x10^'^ grams (where we use a radius of 746 km). However, Castillo- Rogez 
et al. (2007) quote ridge dimensions of 18 km by 200 km over a length of 1600 km, which 
yields a mass of 3x10^^ grams, similar to the Ip (2006) estimate. Profiles from the DTM 
model (Giese et al. 2008; Fig. 5) do not show a recognizable base greater than 50 km for 10 
different profiles across the ridge and no elevations greater than ~13 km, making this a very 
high estimate. 

We assume that the ridge consists of ring material that lands on the surface of lapetus, 
accounting for the equatorial location. Given a ring of material interior to the Roche limit, 
there are two ways for it to land on the surface of lapetus: the ring can tidally evolve down 
to the surface, or it can be pushed there by the newly formed sub-sateUite. For reasons 
discussed below, we focus on the latter. 
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4.2. Sub-satellite and ring interaction 



The tidal evolution of the ring down to the surface of lapetus requires that the material 
be inside the synchronous rotation height. As described above, the synchronous height 
for a rotation period of 16 h is at a = 2.94 Rj, which is exterior to the Roche limit at 
'"roche = 2.53i?i for = 1 g cm~'^. The ring material evolves due to tides with lapetus and 
the sub-satellite that accretes beyond the Roche limit. By comparing our estimates for the 
sub-satellite ( §3.11) to those of the ridge ( §4.11) . we find that the sub-satellite is more massive 
than the ridge for the entire range of sub-satellite masses that assist in de-spinning. In this 
case, the ring spreading timescale i s the time it takes a particle to random walk across a 
distance r iGoldreich and Tremaind (119821 ). 



t. 



spread 



(191 yr) (5600 gem- 



150 kmV 



566 km 



(4) 



S / V -Rss 

where E is the surface density of the ring, and i?ss and a^s are the radius and the semi-major 
axis of the sub-satellite, respectively. The surface density (E) of the ring is simply the ridge 
mass spread over the region interior to the Roche limit, ~ 5.6-46.4 x 10^ g cm~^. The 
possible range of sub-satellite masses is equivalent to the mass of a single body of radius, 
i?ss, of ~131-211 km for a density of 1 g cm~^ (or 155-251 km for P ss=0.6 g cm~^). Th e 



semi-major axis of the sub-satellite is likely to be 1.3 -Rroche initially ( iKokubo et a/.ll2000l ). 
and so the range of possible times for the spreading of the ring is 9-286 years. Thus, even 
with many conservative approximations, this timescale is many orders of magnitude shorter 
than other timescales in the problem. 

For the sub-satellite masses of interest, the effect of the ring on the sub-satellite would 
be dwarfed by the much larger effect of lapetus's tides. 



5. Impact scenarios 

In the vast majority of situations in which the addition of a sub-satellite aids de-spinning, 
the sub-satellite is stripped from its orbit around lapetus. In these cases the stripped sub- 
satellite will still be bound in the Saturnian system, at least initially. In this section, we 
determine the possible fates for these objects and ask what effect this will have on lapetus. 



5.1. The dynamical fate of a stripped sub-satellite 



To determine the probability of impact by a stripped sub-satellite we performed a nu- 
merical A^-body experiment consisting of the orbital evolution of 50 massless test particles 
initially in orbit around l apetus. We used Sy MBA, a symplectic code which is capable of 
handling close encounters iDuncan et al\ (Il998[ ). The simulations included Titan, Hyperion, 
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lapetus, and Phoebe. Particles were stopped if they hit a sateUite, crossed Titan's orbit, 
became unbound from Saturn, or reached a distance of 0.4 AU from Saturn (roughly its Hill 
radius). The particles' initial semi-major axes ranged from 0.4 to 0.8 -Rh (19-38-Ri) with 
eccentricities of 0.1. The particles' inclinations were between 0° and 15° degrees with respect 
to Saturn's equator. 

After 6 Myr only 1 particle remained in orbit around lapetus, while 3 particles hit la- 
petus before becoming unbound. The remaining 46 particles became unbound from lapetus, 
entering orbit about Saturn. These are the objects of interest here because our goal is to de- 
termine the fate of a sub-satellite once it becomes unbound. Five of them were ejected from 
the Saturn system by the satellites (entering heliocentric orbit). The remaining 41 objects 
impacted lapetus — none hit Titan, Hyperion, or Phoebe. Thus, a stripped sub-satellite has 
a roughly 41/46 ~ 90% chance of ending its existence by returning from orbit about Saturn 
and colliding with lapetus. This fact means that we must consider the effects of such an 
impact in our scenario. 



5.2. Angular Momentum Budget 

In cases in which the stripped sub-satellite re-impacts lapetus, we must consider the 
angular momentum it imparts. If the resulting spin rate is too fast it will cancel any ad- 
vantage that was originally gained by the presence of the sub-satellite. We assume that our 
scenario is still viable if, after the impact, fij < fi* = 2 x 10"^, where f2i is the spin frequency 
of lapetus. For spin rates < VL* , corresponding to spin periods greater than ~4 days, the 
de-spinning time after the impact will be less than ~ 100 Myr for ?7=10^^ Pa s. 

Assuming the sub-satellite is accreted completely, the magnitude of the angular momen- 
tum brought in by the sub-satellite is if = nissVoob, where foo is the satellite's velocity with 
respect to lapetus at "infinity", and b is the impact parameter. Assuming that lapetus's 
pre- impact rotation is slow, we find that Qi < Q* requires that b < b* = 2i?jf2*/5gfoo- The 
maximum value of b that allows a coUision with lapetus is 6^^^ = i?j[l + (fesc/^oo)^], where 
fesc is lapetus' surface escape speed. If b* > 6max, all impacts leave lapetus spinning more 
slowly than Q*; otherwise the probability P that lapetus will have Qi < Q* (and thus our 
scenario will remain viable) is (&*/&max)^ = (2i?if2*/5gf )^, where v is the impact speed and 
we have used f ^ = f^^^ + f^. 

In §3.H we found that satellites with 0.005 < q < 0.021 were most effective in de-spinning 
lapetus. Most impacts occurred at velocities near the escape speed of lapetus, 0.58 km/s. 
At that speed, P = 1 for g < 0.0103 and P = (0.0103/^)^ for 0.0103 <q< 0.021. Thus, the 
lower- mass sub-satellites never produce spin rates faster than Q*, while sub-satellites with 
q = 0.0146 yield viable scenarios 50% of the time. We therefore conclude that the re-impact 
of the sub-satellite can be consistent with lapetus's spin state. 
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5.3. Linking basins to the possible sub-satellite impact 

A large complex crater or basin will form if the sub-satellite were to impact lapetus. In 
this section, we estimate the size of the impactors needed to produce the basins observed 
on lapetus today and compare them to the size of the sub-satellite. Zahnle et al. (2003) 
estimate the diameter, D, of the final (collapsed) crater on a mid-sized icy satellite to be 

/ \ 0.49 / 2 \ 0.245 / x 0.377 / „ \ 0.885 

D = 21.4km f-Vl , (5) 



1 km /s/ \ g J \ P J 

where v is the normal component of the impact velocity, g is the surface gravity of the target, 
p is the density of the target, and pi is the density of the impactor. This equation assumes 
that the incidence angle of the impact, measured from the normal to the target, is 45°. If 
we substitute g = 22.3 cm s~^ for lapetus, pi = p = 1.0 g cm~^, we have 

0.49 / p \ 0.885 
^ss 



D = 589 km ^ . (6) 

Vlkm/sy \WOkmJ ^ ' 

Giese et al. (2008) found 7 basins, defined as craters with D > 300 km, on the leading 
face of lapetus. The largest is stated to have D = 800 km. The Gazetteer of Planetary 
Nomenclature (http : / / planetary names . wr . usgs . gov/ ) lists 5 basins on lapetus, with the 
largest, Turgis, 580 km in diameter. 

In our scenario, one of lapetus's basins may have been created by our escaped sub- 
satellite when it returns to meet its maker. In our simulations in §5.H we find that the 
impact would likely occur at a velocity near lapetus's escape velocity, 0.58 km/s; however, 
there is about a 10% chance that the impact speed would exceed 2 km/s. At 2 km/s, Eq. 
indicates that 300 - 800 km basins are produced by impactors with radii between 32 and 
96 km. At the more likely speed of 0.58 km/s, the corresponding impactor radii are 64 and 
193 km. Recall that in §4.2^ we found that our hypothetical sub-satellite should have radii 
between 131 and 211 km for these densities. The ranges vary slightly for a sub-satellite of 
only 0.6 g cm~'^ with impactors between 40-120 km at 2 km/s and 80-241 km at 0.58 km/s, 
where the hypothetical sub-satellite at this density would have a radii between 155-251 km. 
Therefore, it is quite possible that one of the observed basins was caused by our hypothetical 
sub-satellite. 

These calculations are highly uncertain because numerical simulations at the relevant 
scales and velocities have not, to our knowledge, been performed for icy targets. In partic- 
ular, scaling relations in the literature, such as the Schmidt-Housen scaling we used above, 
are generally based on field data or simulations of hypervelocity impacts, i.e., impacts at 
velocities large compared with the speed of sound in the target, which is about 3 km/s for 
non-porous ice. Impacts by lapetus's putative sub-satellite would have occurred at lower 
speeds. Furthermore, experiments and simulations generally deal with impactors that are 
much smaller than their targets. This condition is only marginally satisfied at typical speeds 



- 13 - 



for the impacting sub-satellite in our scenario. Thus, the values we quote for the sub- 
satellite's size should be viewed as rough estimates. However, the calculations do support 
the possibility of basin formation caused by the impacting sub-satellite. 



6. Discussion and conclusions 

We have explored the scenario in which a sub-satellite forms from an impact-generated 
debris disk around lapetus following an impact. The remains of the disk fall to the surface of 
lapetus to build the observed equatorial ridge, while the tidal evolution of the sub-satellite 
assists in de-spinning lapetus. We find that this scenario can significantly shorten lapetus's 
tidal de-spinning time if the mass ratio between the sub-satellite and lapetus, g, is between 
0.005 and 0.021. These results suggest that the presence of a sub-satellite can potentially 
loosen constraints on the geophysical history of lapetus that have been implied by the timing 
and duration of despinning. The full implications of our results, however, cannot be realized 
until the sub-satellite scenario is investigated using a detailed thermal evolution model such 
as those described in prior works (e.g. Castillo- Rogez et al. (2007); Robuchon et al. (2010)). 

lapetus has been of great interest due to its extremely old surface that records the 
cratering history of the outer solar system. Prior scenarios require that the ridge form after 
an epoch of early heating to de-spin lapetus that drive the timing and duration of resurfacing 
(Castillo- Rogez et al. (2007); Robuchon et al. (2010)). In our model, the ridge forms only 
a few hundred years after the impact, and therefore the ridge is quite old. This matches 
the qualitative assessme nt that the ridge is one of the oldest features on the surface, along 



with the 800-km basin ( iGiese et a/.l 120081 ). Some of the sub-satellite evolutionary scenarios 



presented here end with a re- impact, requiring an associated basin forming much later. 

The ridge is only seen to extend roughly 110° around lapetus. This could be the result 
of its extreme age — much of it could have been destroyed by subsequent cratering. Thus, it 
is still unclear whether the ridge extends for the full extent of the equatorial circumference, 
but an infalling ring would preferentially deposit material on global topographic high terrain. 
Indeed, a ring that extends only 110° in longitude could result if the center-qf-figu re of 



lapetus were offset from its center-of-mass as is seen on the Moon (see lAraki et a/.ll2009l . and 
references therein). 

The stripped sub-satellite mass range, described above, corresponds to bodies with radii 
of between roughly 130 and 210 km (assuming pi =1 g cm~^), not all of which are below the 
estimated radius, ~190 km, of the largest allowable impactor given the basins on lapetus. 
This limiting size corresponds to a mass ratio q = 0.015. We arrive at a similar upper limit 
when considering the angular momentum of the impact. Thus, to account for the impact 
of the secondary, the mass ratio range 0.006< q <0.015 is required. As seen in Fig. [21 the 
sub-satellites are stripped on very similar timescales to the de-spinning of lapetus. Thus, a 
sub-satellite returning to form a basin on lapetus would do so after the ridge formed. This 
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younger basin caused by the sub-satellite would then be expected to overlie the ridge if it 
was a near-equatorial impact. Turgis (the 580 km Basin II in Giese et al. 2008) fulfills these 
requirements as the equatorial ridge appears to stop when intersecting this basin (though 
there are no profiles presented in Giese et al. 2008 for this region of the surface). This 
basin would correspond to an impactor with a radius of 140 km, which is only slightly larger 
than the lower limit estimated above. However, we find this acceptable given the inherent 
uncertainties in calculating crater sizes. Meanwhile Basin I, the larges t at 800 km, appears 



to have a stratigraphically similar age as the ridge ( iGiese et a/.ll2008l ). It is important to 
note, however, that our simulations predict only ~90% of stripped sub-satellites return to 
impact, so it is not a certainty that one of the basins is associated with this evolution. 

A strength of our hypothesis is that the limits on impactor and sub-satellite masses 
(0.005 < q < 0.015) are in line with the estimated ridge mass and the number of 300-800 km 
basins. Thus, this work provides a complete story from original to final impact which may 
explain the ridge, shape, and basin population on lapetus. 
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Fig. 1. — The lifetime of each test particle is plotted as a function of their initial semimajor 
axis for two different intial eccentricities (black) 0.1 and (gray) 0.0. The simulations lasted 
for 1 Myr, which is shown as a horizontal line. Symbols for particles which survive for 1 Myr 
are smaller than those of particles with shorter lifetimes. The lifetime drops precipitously at 
a = 21 Ri = 0.44 Rr. 
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Fig. 2. — Results from the integrations of the despinning of lapetus, where despinning 
and sub-satellite stripping/impacting timescales were calculated for different mass ratios q 
for different internal viscosities f] and different internal temperature evolution. The filled 
symbols arc the despinning times for simulations where the sub-satcUitc is stripped, where 
the smaller symbol at the same q is the time of stripping. The open symbols are for cases 
where the sub-sateUite evolves back towards lapetus, and the smaller symbol marks the time 
of re-impact. The times are normalized to the time it takes lapetus to de-spin by Saturnian 
tides alone for each case. Four cases are plotted: (a) constant viscosity, where the results 
from each constant value of 77 = 10^^ Pa s, = 10^^ Pa s, etc. plot on top of each other 
when normalized to 1.0, (b) the LLRI case from Castillo-Rogez (2007), (c) the 0.04 ppb and 
(d) 72 ppb 2^A1 cases from Robuchon et al. (2010). 
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Fig. 3. — Evolution of the spin period of lapetus (solid line) and orbital period of the sub- 
satellite (dashed line) as a function of time, for with constant rj equal to 10^^ Pa s and 
mass ratios of (a) q = 0.024, (b) q = 0.018 and (c) q = 0.004. The dotted line in each 
panel is the evolution of the spin of lapetus under the influence of Saturn only, and the 
shaded region marks the age of the solar system. They illustrate three different outcomes. 
In case (a) the sub-sateUite is too large and eventually is caught in synchronous lock with 
the rotation of lapetus. Saturn continues despinning lapetus, and so the sub-satellite falls 
below synchronous height, returning to impact lapetus. lapetus then despins, again, due to 
Saturn, and finally reaches a despun state later than had it simply despun due to the effects 
of Saturn. In case (b) the sub-satellite assists in despinning lapetus, and is then stripped by 
Saturn, allowing lapetus to dcspin up to 10 x faster than by Saturn alone. Finally, in (c), 
the small sub-satellite's orbit evolves very slowly, so that lapetus is despun by Saturn fast 
enough for the synchronous limit to move beyond the sub-satellite, forcing the sub-satellite 
to tidally contract its orbit and return to impact lapetus. For this case the despinning time 
is similar to that by Saturn's effect alone. 



